Does the isotope effect of mercury support the BCS theory? 
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In this paper, we reexamine the results of isotope effect experiments of the conventional 
monoatomic superconductor (Hg). It is shown clearly that the isotopic coefficients of mercury 
can be largely deviated from a = 0.5, the standard value suggested by the phonon-mediated BCS 
pairing theory. According to the reported experimental results of various mercury isotopes, a giant 
isotope effect (a = 2.896) is numerically found in the data. This study indicates that the validity 
of the conventional BCS theory cannot be verified by the isotope effect. 
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In the field of superconductivity, the isotope effect has 
been considered can play the role in unraveling the mi- 
croscopic superconducting mechanism. Historically, the 
discovery of the isotope effect in mercury in 1950 was 
immediately raised to the question regarding the cou- 
pling of the electrons to the lattice vibrations in the 
superconductor P, It was widely believed that the 
isotope effect may relate to the origin of the effective 
attractive interaction between the repulsive electrons, 
which leads to the occurrence of superconductivity as 
suggested by the BCS pairing theory [3]. In the frame- 
work of BCS pairing theory, it was predicted that there 
would be a universal^ isotopic coefficient a = 0.5 for 
any superconductors [3]. The condensed matter physi- 
cists concluded that this prediction of the BCS theory 
was in good agreement with the reported isotope effects 
in some conventional superconductors (e.g., Hg, Sn and 
Pb)[|. 

The experimental facts show that the isotope effect 
of the vast majority of superconductors (both conven- 
tional and non conventional) largely deviates from the 
standard value 0.5 of BCS theory [6j4l0j. For the conven- 
tional monoatomic superconductors, one notes that some 
superconductors have a negligible coefficient ctzr — 
, while there are some reports on the inverse isotopic 
coefficients 0, for example, <x\j = —2 for the uranium 
element The inverse isotope effect has also been ob- 
served in numerous organic superconductors[10]. For 
the high-temperature superconductors, the isotope ef- 
fect can have coefficients both smaller and larger than 
0.5 depending on the doping levels. One notes that the 
cuprate superconductor Y BaiCu-sO-j-y has a very small 
oxygen isotope effect a — 0.0 ± 0.027 under the tem- 
perature of about 90 Kelvin[9j. Recently, the so-called 
large iron isotope effect was reported in the newly discov- 
ered iron-based SmFeAsOi- x F x and Ba\^ x K x Fe2As2 
superconductors ■ 

In this letter, we will argue for the first time that the 
well-accepted isotopic coefficient of mercury in fact is not 
equal to 1/2. It will be shown that the maximum isotopic 
coefficient of mercury can reach as high as 2.896, a value 



that is about six times of the value predicted by BCS 
theory and widely reported in the related papers. This 
result indicates that the relationship between the criti- 
cal temperature T c the isotopic mass M is much more 
complicated than that of the BCS theory. It seems most 
likely that the value of the isotopic coefficient does not 
directly lead to any priori conclusion about the pairing 
mechanism of the superconductivity. In other words, the 
BCS theory of the electron-phonon interaction mecha- 
nism not only cannot explain the isotope effect of the 
non-conventional superconductors, in fact, it also cannot 
explain the isotope effect observed in the conventional 
systems. 

In the framework of BCS pairing theory, the supercon- 
ducting critical temperature T c is given by 



ksTc = l.lShujoexp 
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VN(Ep) 



(1) 



where cor) is the Debye frequency, V is the elec- 
tron-phonon interaction strength and N(Ef) is the elec- 
tronic density of states at the Fermi surface. 

Eq. ((T|) is considered as one of the most significant 
and influential predictions of the BCS theory. In the har- 
monic approximation, both V and N(Ep) are indepen- 
dent of the ionic mass, while the characteristic frequency 
we can be expressed as 
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U>D cx 



M 



(2) 



where M is the ionic mass. 

By the combination of Eq. ([1} and Eq. ([2]), it is not 
difficult to conclude that if the product VN(Ep) is in- 
creased, then the T c will rise. Moreover, using light ele- 
ments M can raise the Debye frequency , and this in 
turn enhances the T c of the corresponding superconduc- 
tors. Based on the above two equations, BCS predicted 
a maximum T c of around 30K for any superconductors 
which has been proven wrong since the discovery of the 
high-T c oxide superconductors by Bednorz and Miiler in 
19860. 
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Figure 1: The experimental results of superconducting tran- 
sition temperature as a function of isotopic mass for the ele- 
mental mercury. 



Within the framework of the electron-phonon mecha- 
nism, the T c can be described by the following relation: 



T c = AM~ 



(3) 



where A is a constant, M is the mass of the element 
substituted by its isotope and a is the so-called isotope 
coefficient which is defined as 



d\nT c 
"ainM 



M AT C 
"f, AM' 



(4) 



where AT C is the shift of the critical temperature and 
AM is the difference between the two isotopic mass. 

In the standard BCS theory, T c is inversely propor- 
tional to the square root of the masses of the isotopic el- 
ements, hence the isotope-effect coefficient a = 0.5 which 
has been considered in good agreement with the experi- 
mental results in some non-transition metal superconduc- 
tors, such as Hg, Sn and Pb. In the following discussion, 
we will show that this well-known conclusion is evidently 
wrong. 

The superconductivity was first observed in mercury 
by Onnes in 191l[l2j. Later, the isotope effects for mer- 
cury had been intensively investigated and its isotope- 
effect was claimed to be around 1/2, as shown in Fig. 
[TJ These results were the basis and foundation of the 
electron-phonon mechanism of BCS theory. And now, 
in order to illustrate the validity of the BCS theory, the 



Table I: The reported superconducting transition tempera- 
ture T c (i) of mercury with an isotopic mass of Mj. Moreover, 
at and /3i are newly defined parameters of isotope-effect coef- 
ficient and the percent error, respectively. 



Hg(i) 


1 


2 


3 


4 


5 


6 


Mi 


198.0 


199.7 


200.6 


200.7 


202.0 


203.4 


T c (i) 


4.177 


4.161 


4.156 


4.150 


4.143 


4.126 




0.446 


0.267 


2.896 


0.260 


0.592 




Pi 


10.8 


46.6 


479.2 


48.0 


18.4 





figure has been widely adopted in the textbooks of su- 
perconductivity. 

Now let us begin our study of the relationship between 
the isotope effect and the BCS theory. For the sake 
of convenience of discussion in this letter, the isotopic 
masses Mi and the corresponding superconducting tran- 
sition temperatures T c {i) are listed in the table below. 
It became obvious that only for the purpose of demon- 
strating the correctness of the prediction a = 0.5 of BCS 
theory, the researchers intentionally selected special iso- 
topic samples (1 and 6) and the isotope-effect coefficient 
was calculated by 



Mi 



T c (6) - T c (l) 



M 6 - Mi 



0.448. 



(3) 



With the two sets of data (T c (l), M x ) and (T c (6), M 6 ), 
one indeed can obtain the desired isotope value close to 
1/2, but this is completely wrong. Mathematically, to 
apply Eq. (j4|) to estimate the parameter a, the variable 
quantities AT C and AM must be infinitesimal. Hence, 
the calculation formula ((5|) holds only when AT C and 
AM satisfy a linear relationship which is obviously not 
consistent with Eq. ([3]). Here, it is argued that a more 
accurate formula for calculating the isotope-effect coeffi- 
cient would be 



lnT c (i + l)-lnr e (z) 

In Mi+i - In Mi 
Mj T c (i + 1) - T c (i) 
T c (i) M i+1 -Mi 



(G) 



In the above formula, two sets of adjacent data (T c (i), 
Mi) and (T c (i + 1), Mj +1 ) are applied in the numerical 
simulation. We have argued that for a given supercon- 
ductor, its isotope-effect coefficient may be completely 
different. This argument is well confirmed by directly 
putting the reported experimental results of Fig. [T] into 
Eq. IS]) , as shown in Table |TJ It is not difficult to find 
that all the obtained m are very different from 0.5 (in the 
range from 0.2 to 3), surprisingly, a giant isotope effect 
(a = 2.896 ) which is about six times higher than the 
value 1/2 suggested by BCS theory. In order to better 
present the deviation of the experimental results from the 
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BCS theory, we define the percent error ft as: 

ft = x 100%. (7) 

where 0.5 is the full isotope effect in the framework of 
BCS theory. 

As also shown in TableU the minimum percent error of 
the isotope coefficient for the elemental mercury is about 
10.8%, while the maximum percent error can reach an 
inconceivable value of 479.2%. These results imply that 
the validity of the full isotope effect (a = 0.5) of BCS 
pairing theory has never been experimentally verified. It 
is a possibility that the electron-phonon interaction based 
BCS theory may be fundamentally flawed. 

In this short letter, we have presented arguments 
against the mainstream view that the isotope effect in 
the non-transition metal superconductor (Hg) is equal to 
1/2, which is the characteristic value predicted by the 
classical form of BCS theory. It has been pointed out 
that the traditional method of calculation used to es- 
timate superconducting isotope effect is mathematically 
unreliable. Based on the reported experimental result, 
very different isotope coefficients ranged from 0.267 to 
2.896 for the mercury have been numerically obtained 
by using the improved method. Our results imply that 
the dependence of the critical temperature for supercon- 
ductivity upon the isotopic mass was much more compli- 
cated than the BCS theory previously advocated. In our 



opinion, the isotope effect cannot be applied as the di- 
rect evidence of the proposed electron-phonon coupling 
mechanism. 
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